A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

02.05.2; 02.10.8; 08.10.1; 08.14.1; 08.16.6 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OO 
ON 

On 
PL. 



> 
(N 

m 

o : 

00 . 
On ■ 

^ ' 

Or 

6 
c3 



Equilibrium Sequences of Rotating Neutron Stars for New 
Realistic Equations of State 

B. Datta 1 0, A. V. Thampan^nd I. Bombaci 2 

1 Indian Institute of Astrophysics, Bangalore 560 034, INDIA. 

2 Dipartimento di Fisica, Universita di Pisa, and INFN Sezione di Pisa, Piazza Torricelli 2, 1-56100 Pisa, ITALY 



Abstract. For four newly suggested realistic equations of state 
of neutron star matter, we construct equilibrium sequences 
of rapidly rotating neutron stars in general relativity. The 
sequences are the normal and supramassive evolutionary se- 
quences of constant rest mass. We find that for these equa- 
tions of state the maximum (gravitational) mass rotating mod- 
els occur (in central density and rotation rate Q) before the 
maximum- f2 models. We calculate equilibrium sequences for 
a constant value of Q corresponding to the most rapidly ro- 
tating pulsar PSR 1937+21. Also calculated is the radius of 
the marginally stable orbit and its dependence on Q, relevant 
for modeling of kilo-Hertz quasi-periodic oscillations in X-ray 
binaries. 
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1. Introduction 

Equilibrium sequences of rapidly rotating neutron stars are im- 
portant in modelling a variety of phenomena of astrophysical 
interest, such as millisecond pulsars, low-mass X-ray binaries 
(LMXBs) and Quasi Periodic Oscillators (QPOs). Models of 
rapidly rotating neutron stars in general relativity must be 
constructed numerically. Early work on this have been based 
on incompressible fluids and polytropic models (Bonazzola & 
Schneider 1974; Butterworth 1976). In 1986 Friedman et al. 
(1986) reported calculations of rapidly rotating neutron stars 
in general relativity using a set of realistic equations of state 
(EOS) for neutron star matter. A similar work based on a for- 
malism due to Komatsu et al. (1989) (KEH formalism) was 
done by Cook et al. (1994) for purpose of studying quasi- 
stationary evolution of isolated neutron stars. An alternative 
approach based on spectral methods was developed by Bonaz- 
zola et al. (1993). Extensive calculations using the spectral 
method for a broad set of realistic EOS of neutron star matter 
were presented in Salgado et al. (1994a; 1994b). 

A key input in determining the structure of neutron stars 
is the EOS of high density matter. The work of Friedman et al. 
(1986) and that of Cook et al. (1994) make abundantly clear 
that the EOS also plays a significant role in deciding the various 
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equilibrium sequences of rotating neutron stars. For example, 
the Keplerian frequency of a test particle in orbit around a 
neutron star ranges from 55% of its spherical value for models 
based on the softest EOS to 75% of the spherical value for 
models with the stiffest EOS. The spreads in rotation-induced 
changes in the values of masses and radii from static neutron 
star cases also display considerable EOS dependence. These 
quantities (especially the Keplerian frequency of a particle in 
orbit around the rotating neutron star and the radius of the 
innermost stable circular orbit) are important for deciding the 
boundary layer structure and hence the emission characteristic 
of LMXB and QPOs. 

Although dense matter has been a subject of study for 
nearly three decades now, there is no general agreement still 
on its exact composition, and on its EOS, especially for den- 
sities in excess several times nuclear matter density. The bulk 
of a neutron star (core) is made up of an electrically neutral 
quantum fluid composed of neutrons, protons, electrons and 
muons in equilibrium with respect to the weak interactions 
(beta-stable nuclear matter). However, at ultra-high densities, 
a variety of new and exotic hadronic degrees of freedom may 
become important (like hyperons, a K~ condensate or a decon- 
fined phase of quark matter). The possible appearance of such 
an exotic core has enormous consequences for the transport 
properties of neutron stars and also for the formation of black 
holes (Brown & Bethe 1994; Bombaci 1996). The consequences 
of the existence of an exotic core (such as quark matter or 
kaon condensation) on the properties of rapidly rotating neu- 
tron stars will be reported in a forthcoming paper (Datta et al. 
1997). In the present work, we have considered a conventional 
picture assuming the neutron star core to be composed of only 
beta-stable nuclear matter. Even in this picture, the determi- 
nation of the EOS of asymmetric nuclear matter to describe 
the core of the neutron star, remains a formidable theoretical 
problem. In fact, one has to extrapolate the EOS to extreme 
conditions of high density and high neutron-proton asymme- 
try, i.e. in a regime where the EOS is poorly constrained by nu- 
clear data and experiments. Astrophysical observational data, 
such as based on the binary pulsar PSR 1913+16, which give 
1.4 Mq as the mass of the neutron star in the binary (Taylor 
& Weisberg 1989) and analysis of Vela pulsar postglitch tim- 
ing data can be used to broadly rule out very soft EOS (Datta 
& Alpar 1993). Recently, some new EOS of asymmetric nu- 
clear matter have been calculated and applied to the study of 
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non-rotating neutron stars (Baldo et al. 1997, Bombaci 1995). 
These EOS are based on (i) a microscopic Brueckner-Bethe- 
Goldstone many-body approach and (ii) a phenomenological 
model based on effective nuclear forces. These satisfy the basic 
requirements of reproducing the empirical saturation point for 
symmetric nuclear matter, the symmetry energy and the in- 
compressibility parameter at the saturation density (see Table 
1). These models have the desirable physical feature that the 
velocity of sound in the medium does not violate the causality 
condition. Therefore, these can be taken to be realistic EOS, 
and so it would be of interest to see the equilibrium rotating 
sequences that would be possible with these EOS. In this paper 
we report calculations of equilibrium sequences of rapidly ro- 
tating neutron stars in general relativity for these new realistic 
EOS models. The various equilibrium sequences that we con- 
struct are normal and supramassive evolutionary sequences of 
constant rest mass. In addition, we build equilibrium sequences 
for a constant value of rotation rate corresponding to a pe- 
riod of P = 1.558 ms of the millisecond pulsar PSR 1937+21 
(Backer et al. 1982), the most rapidly rotating pulsar known. 

2. Rapidly and Rigidly Rotating Relativistic Stars 

The space-time around a rotating neutron star can be de- 
scribed in quasi-isotropic coordinates, as a generalization of 
Bardeen's metric (Bardeen 1970): 

ds 2 — g fJ , l/ dx p dx v (p, v — 0, 1, 2, 3) 

= -e /+p dt 2 + e 2a (r 2 d6 2 + dr 2 ) + e 1 - p r 2 sin 2 9 

(d<p-udtf (1) 

where g M „ is the metric tensor. The metric potentials 7, p, a, 
and the angular velocity of the stellar fluid relative to the local 
inertial frame (u) are all functions of the quasi-isotropic radial 
coordinate (r) and the polar angle (9). We use here geometric 
units: c = 1 = G. Wc assume a perfect fluid description, for 
which the energy momentum tensor is given by: 



Substitution of the above into Einstein field equations pro- 
jected on to the frame of reference of a ZAMO yield three 
elliptic equations for the metric potentials p, 7 and ui and two 
linear ordinary differential equations for the metric potential 
a (Komatsu et al. 1989; Butterworth & Ipser 1976; Bardeen & 
Wagoner 1971). In the KEH formalism (Komatsu et al. 1989), 
the elliptic differential equations are converted to integral equa- 
tions for the metric potentials using Green's function approach. 

From the relativistic equations of motion, the equation of 
hydrostatic equilibrium for a barytropic fluid may be obtained 
as: 

h(P) -h p = / ^- = km* - km* - / F(Q)dil (6) 

where h(P) is termed as the specific enthalpy. P p , u* are the 
rescaled values of pressure and t-component of the four velocity 
respectively and h p is the specific enthalpy at the pole; F(Q) = 
u Utj, is the integrability condition imposed on the equation of 
hydrostatic equilibrium, and it can be physically interpreted 
as the rotation law for the matter constituting the neutron 
star. An appropriately chosen value of hp defines the surface of 
the star. Equation (6) shows that the hydrostatic equilibrium 
equation is integrable if P(e) and u are specified. 

As shown by Bardeen 1970 (see also Butterworth & Ipser 
1976), the quantity is a function of f2 only. Komatsu et 

al. (1989) have suggested the following specific form for F(Q): 

F(Q) = A 2 (Q c -n) (7) 

where A is a rotation constant such that when A — > 00, the 
configuration approaches rigid rotation (that is, — Q c ) so 
that F(Q) is finite. Furthermore, when A — > 0, the configu- 
ration should approach that of rotation with constant specific 
angular momentum. 

On substituting Eqs. (5) and (7) into Eq. (8), we have the 
hydrostatic equilibrium equation as 



= (e + P)// + Pg pv (2) 

where e is the total energy density, P the pressure and u p the 
unit time-like four velocity vector that satisfies 

/m, = -1 (3) 

The proper velocity v of the matter, relative to the lo- 
cal Zero Angular Momentum Observer (ZAMO), is given in 
terms of the the angular velocity Q, = u 3 /u° of the fluid ele- 
ment (measured by a distant observer in an asymptotically flat 
space-time), by the following equation (see Bardeen 1970): 

v = (n~Lo)rsinde- p (4) 
The four velocity (u M ) of the matter can be written as 

-( 7 +p)/2 

n-i.W 1 - - -") (5) 



h(P) -h p = - [7 P + p p - 7 - p - ln(l - v 2 ) + 

A 2 (n-n c ) 2 ] (8) 

where 7 P and p p are the values of the metric potentials at the 
pole, and Q — r e Q.. 

Therefore, the hydrostatic equilibrium equations at the 
centre and equator for a rigidly rotating neutron star become 
respectively 

h(P c ) - hp - i [7 P +p p - 7c -p c ] = (9) 

(7p + P P -le- p e ) - ki[l - (n e - cu e ) 2 r 2 e- 2p "] = (10) 

where the subscripts p, e and c on the variables stand respec- 
tively for the corresponding values at the pole, equator and 
center. 

We solve (numerically) the integral equations for p, 7 and 
u), the ordinary differential equation (in 9) for the metric po- 
tential a, together with Eqs. (8), (9) and (10), iteratively to 
obtain p, 7, a, w, the equatorial coordinate radius (r e ), angular 
velocity (O), and the density (e) and pressure (P) profiles. 
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2.1. Innermost stable orbits 

Since the metric (1) is stationary and axisymmetric, the energy 
and angular momentum are constants of motion. Therefore, for 
a particle in stable orbit around the neutron star, the specific 
energy E (in units of the rest energy moc 2 , where mo is the 
rest mass of the particle) and the specific angular momentum I 
(in units of moc) can be identified as — p and p 3 respectively, 
where, p M (/x = 0, 1, 2, 3), stands for the four-momentum of the 
particle. From the condition p M p M = — 1, we have the equations 
of motion of the particle (confined to the equatorial plane) in 
this gravitational field as 



; dt o -H+P) 
t = — — p = e 1 ' H> 
dr 



(E-u>l) 



•2 ^ e 2 Q+7+P ^ = E 2_y2 



(11) 

(12) 
(13) 



Here, r is the proper time and V is the effective potential given 
by 



V 2 = e 1+p 



1 + 



l 2 /r 2 

e i-p 



+ 2coEl~u> 2 l 2 . 



(14) 



The conditions for circular orbits, extremum of energy and 
minimum of energy are respectively: 



E 2 = V 2 

V, r = 
V rr > 0. 



For marginally stable orbits, 



V, r 



0. 



(15) 
(16) 
(17) 



(18) 



In our notation, a comma followed by one V represents a first 
order partial derivative with respect to r, etc.. 

From the expression for the effective potential and the con- 
ditions (15), (16) and (17), one obtains three equations in the 
three unknowns: namely, r, E, and I. In principle, if analytical 
expressions for e 7+p , e 2a , e 7 ~ p and ui are known, it would be 
a straightforward exercise to solve these equations to obtain 
r, E, and I. In practice, however, this is not so, and the solu- 
tions for the metric coefficients e 7+p , e 2a , e 7 ~ p , and ui have to 
be obtained as arrays of numbers for various values of r and 
6 using a numerical method. Furthermore, the condition (18) 
will introduce second order derivatives of 7, p, and lo, which 
means that care has to be exercised in ensuring the numerical 
accuracies of the quantities calculated. For this purpose, it is 
convenient to express E and I in terms of the physical velocity 
v using Eq. (4) (Bardeen 1972) as: 



E-ujI = 



e (7+p)/2 

Vl-v 2 



(19) 
(20) 



Eqs. (19) and (20) can be recognized as the condition for cir- 
cular orbits. Conditions (16) and (18) yield respectively, 



(^/e-^r-Vr + 2r( 7 , r + p, r ) + r 2 ( 7 2 r - p 2 r )z 



p r 2 LO >r ^j/(2 + r(j, r - P,r) 



r I 2 



2 \ 
7,r) - 



1 -1p 2 3 

-e uj r r - 



P,r + 



+ [2 + r(7, r — p,r)]vv ; r — e p uj. r rv 



+ 2 <CI ' r ~ '''' ' 



e p r 2 uj, r v,r = 



(21) 



(22) 



where we have made use of Eq. (21) and its derivative with 
respect to r in order to eliminate the second order derivatives 
in Eq. (22). The zero of V, rr will give the innermost stable cir- 
cular orbit radius (r or b). In Eq. (21), the positive sign refers to 
the co-rotating particles and the negative sign to the counter- 
rotating particles. In this study we have considered only the 
co-rotation case. 

3. Numerical Procedure 

The numerical procedure followed by us is the KEH formal- 
ism. This is based on an earlier work by Hachisu (1986) which 
has a self-consistency requirement that requires that the max- 
imum (central) energy density e c and the ratio of the polar to 
equatorial radial coordinates r p /r e be fixed for each iterative 
cycle. If p\ Y , a 1 and r* are the values of the corresponding 
parameters during the i th iterative cycle, then: 

1. these values are first scaled (divided) by (r' e ) 2 to obtain 
p l , Y and a 1 respectively. 

2. a new value of r e is calculated using Eq. (10) for e — e c 
i.e. v — so that 



2[h(P(e c )) - h, 



the value of S7 C is computed from Eq. (11) as 



Q 



i + l 



u e + e p 



(23) 



(24) 



the values of the three scaled metric potentials p 1 , 7 1 
and a 1 are rescaled (multiplied) by (r* +1 ) 2 
using these values of r* +1 , , p l , 7 \ a 1 , Cj % , equation 
(9) is solved to obtain the matter energy distribution 
namely P l+1 , v t+1 etc. 

the integral equations for the metric potentials are 



solved to obtain p 



+1 



r 



and a l+1 . 



7. steps (1) to (6) are repeated until r e converges to within 
a tolerance of 10 -5 . 

Once r e converges, the metric potentials p, 7, ui and a together 
with the density (e) and pressure (P) profiles can be used to 
compute the structure parameters (see Cook et al. 1994). 

4. New Equation of State Models 

4-1. Microscopic equation of state 

In a microscopic approach the input is the two-body nucleon- 
nucleon (NN) interaction, described by so-called realistic in- 
teractions like the Argonne, Bonn, Nijmegen, Paris, Urbana 



4 



potentials (see e.g. Machleidt 1989). The theoretical basis to 
construct these realistic NN potentials is the meson-exchange 
theory of nuclear forces. In this scheme, nucleons, nucleon res- 
onances (e.g. A(1232)), and mesons such as 7r, p and w, are 
incorporated in a potential representation. The various param- 
eters in the potential are then adjusted to reproduce the exper- 
imental data for the two-body problem (deuteron properties 
and NN scattering phase shifts). Then one has to solve the 
complicated many-body problem to get the EOS. 

Recently, Baldo et al. (1997), hereafter BBB, have com- 
puted a new EOS of beta-stable nuclear matter, and with this 
EOS they have calculated the structure of non-rotating neu- 
tron stars. In their approach, the energy per nucleon of nuclear 
matter is obtained in the Brueckner-Hartree-Fock (BHF) ap- 
proximation of the Brucckncr-Bcthe-Goldstone theory. The 
only input quantity for these calculations is the nuclear inter- 
action. In their calculation BBB used the Argonne vl4 (Avl4)) 
(Wiringa et al. 1984) or the Paris (Lacombe et al. 1980) two- 
body nuclear force, implemented in both cases by the Urbana 
three-body force (TBF) (Carlson et al. 1983; Schiavilla et al. 
1986). As is well known, the need for a TBF arises to ob- 
tain a correct saturation point of symmetric nuclear matter 
in a non-relativistic many-body approach. In the following 
we refer to the EOS obtained in BBB with the Avl4+TBF 
and Paris+TBF, as the BBB1 and BBB2 equation of state re- 
spectively. The saturation properties for these two microscopic 
models are summarized in Table 1, and the calculated speed 
of sound is shown in Fig. 1. The latter always remains within 
the causality bound. 

4-2. Phenomenological equation of state 

In this case the input is a density-dependent effective NN in- 
teraction. The most popular of this kind of interaction is the 
Skyrme interaction (Skyrme 1956). In the present work we 
used a generalized Skyrme-like EOS developed in ref. (Bom- 
baci 1995), and we refer to it as the BPAL EOS. An important 
feature of the BPAL models is the possibility to have different 
forms for the density dependence of the potential part E^y^n) 
of the nuclear symmetry energy, modelling different results pre- 
dicted by microscopic calculations (Wiringa et al. 1988, Bom- 
baci & Lombardo 1991). In particular, E^"^ is proportional 
to the nucleon number density n in the case of BPAL32, and 
to ^Jn in the case of BPAL21 EOS. The density dependence of 
the symmetry energy plays a very important role in the physics 
of neutron stars. This function determines the proton fraction 
in beta-stable nuclear matter, which, in turn, is crucial for an 
accelerated rate of cooling of a neutron star through the so- 
called direct Urea process above a critical value of the proton 
fraction (Lattimer et al. 1991; Page & Applegate 1992). 

In Fig. 2 we plot the pressure-density relationship for EOS 
models BBB1, BBB2, BPAL21 and BPAL32 (curves 1-4) and 
also compare it with two other EOS models, one of which is a 
very soft EOS BPAL12 (Bombaci 1995) and the other a very 
stiff EOS (Sahu et al. 1993): 

(1) Soft: BPAL21 {Bombaci 1995): This EOS is character- 
ized by K = 120 MeV and E^* m ~ n. The value 120 MeV 
for the incompressibility is unrealistically small when com- 
pared with the value 220 ± 30 MeV extracted from nuclear 
phenomenology (Blaizot 1980; Myers & Swiatecky 1996), how- 
ever BPAL12 EOS is still able to sustain the measured mass 
1.44 Mq of the pulsar PSR1916+13 as the maximum gravi- 



tational mass of non-rotating neutron stars constructed with 
this EOS is 1.46 M . 

(2) Stiff: Sahu et al. (1993): This is a field theoretical EOS 
for neutron matter in beta equilibrium based on the chiral 
sigma model. The model includes an isoscalar vector field gen- 
erated dynamically and reproduces the empirical values of the 
nuclear matter saturation density and binding energy and also 
the isospin symmetry coefficient for asymmetric nuclear mat- 
ter. The energy per nucleon of nuclear matter according to 
Sahu et al. (1993) is in very good agreement, up to about four 
times the equilibrium nuclear matter density, with estimates 
inferred from heavy-ion collision experimental data. The max- 
imum gravitational mass of non-rotating neutron stars con- 
structed with this EOS is 2.59 M© . 

For our computations, we constructed the composite EOS 
for the entire span of neutron star densities by joining the new 
high density EOS models to that of Negele & Vautherin (1973) 
for the density range (10 14 — 5 x 10 10 ) g cm~ 3 , Baym et al. 
(1972) for densities down to ~ 10 3 g cm -3 and Feynman et al. 
(1949) for densities less than 10 3 g cm~ 3 . 

5. Results and Discussion 

The equilibrium sequences of rotating neutron stars depend 
on two parameters: the central density (e c ) and the rotation 
rate (fl). For purpose of illustration, we choose three limits in 
this parameter space. These are: (i) the static or non-rotating 
limit, (ii) the limit at which instability to quasi-radial mode 
sets in and (iii) the centrifugal mass shed limit. The last limit 
corresponds to the maximum f2 for which centrifugal forces are 
able to balance the inward gravitational force. 

Table 2 summarizes the non-rotating neutron star struc- 
ture parametes for the EOS models BBB1, BBB2, BPAL21 and 
BPAL32. The values listed correspond to the maximum stable 
mass configuration. The entries in this table are the central 
density (e c ), the gravitational mass (Mg), the rest (baryonic) 
mass (Mo) of the neutron star and the radius (R). The maxi- 
mum mass is an indicator of the softness/stiffness of the EOS 
and its values as listed in Table 2 reflect that the EOS models 
used are all intermediate in stiffness. Among these, BPAL21 is 
the softest EOS and BPAL32 the stiffest one. 

In Table 3 we list the following quantities corresponding to 
the maximum gravitational mass configurations: central den- 
sity, rotation rate, moment of inertia (7), gravitational mass, 
ratio of rotational kinetic energy to total graviational energy 
(T/W), equatorial radius, eccentricity (e) where e is defined 
to be y/\ — Rp/R 2 , R p being the polar radius of the config- 
uration, ratio of rotation rate (u> c ) of stellar fluid relative to 
the inertial frame at the centre of the star to the rotation rate, 
the angular momentum (J), the value of the radius of the in- 
nermost stable orbit (r or b), the polar, forward and backward 
redshifts (Z p , Zf, Zj,), the rest mass (Mo) and the proper mass 
(M p ) of the neutron star. The values of these quantities are 
listed in Table 4 for the maximum angular momentum models. 
From Table 3, it can be seen that the gravitational mass of the 
maximum stable rotating configuration has a value that is close 
to 2 Mq . Interestingly, this is close to predictions from analysis 
of LMXB observational data (Zhang et al. 1996). Therefore, if 
the internal constitution of the compact star in LMXBs were to 
be described by the EOS models that we have considered here, 
and furthermore, if it were to have a mass ~ 2 Mq , the star 
has to be rotating near centrifugal break-up speeds (rotation 
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periods ~ 0.5 ms). Interestingly, for such configurations, it can 
be seen from Table 3, that the separation of the innermost sta- 
ble orbit from the neutron star surface (namely, the boundary 
layer extent) is highly EOS dependent. This separation ranges 
from 0.07 km to 1.031 km. This relatively large spread can 
be understood in terms of the spreads on values of e c and SI 
corresponding to the centrifugal mass shed. These results will 
have relevance in modeling LMXBs/QPOs. For such neutron 
stars as the central objects in LMXBs, there will be a very 
significant re-ordering of the contributions of the disk and the 
boundary layer luminosities to the total luminosity (Thampan 
& Datta 1997). 

The results of our computations for rotating neutron stars 
corresponding to the four new EOS models are given below. 



5.1. EOS Model BPAL21 

The normal sequences for this EOS have rest mass Mq/Mq < 
1.9395 and supramassive sequences have rest mass 1.9395 < 
Mo/Mq < 2.2515. 



5.2. EOS Model BBB1 

In Fig. 3 we show the functional dependence of the gravita- 
tional mass with central density. In this and all subsequent fig- 
ures, the bold solid curve represents the non-rotating or static 
limit, and the bold dashed curve the centrifugal mass shed 
limit. The long dashed curve is the constant-f2 sequence cor- 
responding to the period P — 1.558 ms. The thin solid curves 
that are roughly horizontal are the constant rest mass evolu- 
tionary sequences. The evolutionary sequences above the maxi- 
mum stable non-rotating mass configuration are the supramas- 
sive evolutionary sequences, and those that are below this limit 
are the normal evolutionary sequences. The almost vertical 
thin dashed line is the limit for instability against quasi-radial 
modes. The supramassive evolutionary sequences beyond the 
quasi-radial mode instability limit are represented by dotted 
lines. The maximum mass sequence for this EOS corresponds 
to a rest mass value of 2.471 Mq . The supramassive sequences 
lie in the rest mass range of 2.356 Mq < M < 2.471 M© . 
The gravitational mass of the maximum stable rotating config- 
uration is 2.135 Mq and its radius is 13.129 km. If we assume 
that the fastest pulsar known to date, PSR 1937+21, has the 
canonical mass value of 1.4 Mq and is described by EOS model 
BPAL21, then this neutron star should have a central density 
of about 1.2 x 10 15 g cm' 3 . 

In Fig. 4 we give a plot of Mq as a function of R. For the 
millisecond pulsar PSR 1937 + 21 with an assumed mass value 
of 1.4 Mq , this corresponds to a radius of 11 km. 

In Fig. 5 we display the plot of SI as a function of the spe- 
cific angular momentum cJ/GMq . The inset shows a close-up 
view of the region surrounding the instability limit to quasi- 
radial mode near the centrifugal mass shed limit. It is clear 
from this figure that the maximum mass rotating model (rep- 
resented by the plus sign) has a lower angular velocity than the 
maximum-^ model (represented by the intersection of the line 
representing the instability to quasi-radial modes with that of 
the centrifugal mass shed limit). 



5.3. EOS Model BBB2 

The equilibrim sequences for BBB2 are displayed on Figs. 6 
- 8 for the same representative set of parameters as for EOS 
model BBB1 

For this EOS, the supramassive sequences have rest masses 
between 2.261 Mq and 2.653 Mq and the maximum mass 
at mass shed limit is 2.272 Mq with an equatorial radius of 
12.519 km. 

5.4. EOS Model BPAL32 

This EOS model being the stiffest out of the four models that 
we consider here, has the highest value for the maximum ro- 
tating gravitational mass (2.3 Mq ). The rest masses of supra- 
massive sequences lie in the range 2.263 Mq < Mo < 2.655 

Mq . 

Since the behaviour of M with e c and R and that of SI with 
cJ/GM' 2 for EOS models BPAL21 and BPAL32 are more or 
less similar to those for EOS models BBB1 and BBB2, we 
do not display the corresponding figures for the former EOS 
models here. 

In Table 5 we list the values of the various parameters for 
the constant SI sequences for the four EOS models considered 
in this work. In general, r orb exhibits three characteristics: (a) 
r or b is non-existent (b) r orb < R, and (c) r or b > R. For the 
first two cases, r or t is taken to be the Keplerian orbit radius 
at the surface of the star. From Table 5 it can be seen that 
for low central densities, stable orbits can exist all the way up 
to the surface of the neutron star but for high enough central 
densities, the boundary layer (the separation between the sur- 
face of the neutron star and its innermost stable orbit) can be 
substantial (~ 5 km for the maximum value of the listed cen- 
tral densities) . These results will have applications in modeling 
accretion flows in LMXBs. 

We now make a brief reference to other similar work. For 
equilibrium Keplerian angular velocity corresponding to the 
period of millisecond pulsar PSR 1937+21 and an assumed 
mass of 1.4 Mq for the neutron star, Friedman et al. (1986) 
suggest that stiff EOS for neutron star interior are favoured. 
A similar conclusion but based on a pulsar glitch model and 
the crustal moment of inertia considerations has been reported 
by Datta & Alpar (1993). The work of Friedman et al. (1986) 
show that for a given EOS, the models with maximum grav- 
itational mass also have the greatest frequency of rotation. 
Cook et al. (1994) found that while models with maximum 
gravitational mass also (due to stability conditions defined by 
Friedman et al. 1988) have the maximum rotation rate SI, the 
models for maximum gravitational mass and maximum-^ do 
not in general coincide. In particular, for EOS models that 
display causality violation near or before the maximum sta- 
ble mass non-rotating configuration, the maximum-fi model 
occurs before (in central density and SI) the maximum mass 
model at the mass shed limit. The EOS models that we have 
considered here do not violate the causality condition until 
well beyond the maximum stable mass non-rotating configura- 
tion. Our computations show that the maximum gravitational 
mass rotating models for these EOS occur (in central density 
and SI) before the maximum-fi models. In view of absence 
of correlation between QPO frequency and source count rate 
(as suggested by recent observations, see Berger et al. 1996; 
Zhang et al 1996), our results on marginally stable Keplerian 
orbits corresponding to realistic EOS will have application in 
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understanding kilo-Hertz quasi-periodic oscillations in X-ray 
binaries in terms of strong-field general relativity, where rapid 
rotation of the accreting neutron star is important. 



Zhang W., Lapidus I., White N. E., Titarchuk L., 1996, ApJ 
469, L17 
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TABLE CAPTIONS 

TABLE 1: Saturation properties for the various equations of state used in the 
present work: n [fm~ 3 ] is the saturation number density, E / A[M eV / fm 3 ] the 
corresponding energy per nucleon and A" [MeV] the incompressibility of sym- 
metric nuclear matter. E sym (n )[M eV] is the symmetry energy at the saturation 
point. In the last row we report the empirical saturation properties of nuclear 
matter. 

TABLE 2: Maximum mass non-rotating models: the listed quantities are EOS 
models, central density (e c ) in units of g cm~ 3 , gravitational mass (Ma) in 
solar units, radius (R) in km, baryonic mass (M ) and proper mass (Mp) in 
solar units and lastly the radius (r or b) of the innermost stable orbit in km. Note: 
the numbers following the letter E in column 2 of this and the subsequent tables, 
stand for powers of ten. 

TABLE 3: Maximum mass rotating models: the listed quantities are EOS mod- 
els, central density (e c ) in units of g cm~ 3 , angular velocity (0) in units of 
10 4 rad s~ x of the neutron star as measured by an observer at infinity, moment 
of inertia (/) of neutron star in units of 10 45 g cm 2 , gravitational mass (Mq) 
in solar units, ratio of the rotational energy to the total gravitational energy 
(T/W), radius (R) in km, eccentricity (e), ratio of inertial frame dragging at 
the center of the star to the rotation rate (cu c /O c ), angular momentum (J) in 
10 49 cgs units, radius (r or b) of the innermost stable orbit in km, the polar (Z p ), 
forward (Zj) and backward (Z&) redshifts, baryonic mass (M ) and proper mass 
(Mp) in solar units. 



TABLE 4: Maximum angular momentum models 

TABLE 5: Constant sequences: the sequences have angular velocity of the fastest 
known pulsar PSR 1937+21 



FIGURE CAPTIONS 

FIG 1: Speed of sound in beta-stable nuclear matter as a function of the nucleon 
number density for different models of EOS used in the present work. The three 
continuous lines refer from top to bottom of the figure to the BBBI, BPAL32 
and BPAL2I model respectively; the dashed line is relative to the BBB2 EOS. 
The dot on each curve gives the speed of sound at the center of the non-rotating 
maximum mass configuration for that EOS model. 

FIG 2: Pressure as a function of density for the new EOS models. Curves 1-4 
stand for BBBI, BBB2, BPAL21 and BPAL32 respectively. Also shown for 
comparison are the EOS models BPAL12 (Bombaci 1995) a very soft EOS (rep- 
resented by BPAL12 in figure) and Sahu, Basu & Datta (1993) a very stiff EOS 
(represented by CSM in figure). 

FIG 3: The functional dependence of the gravitational mass with central density 
for EOS model BBBI. In this and all subsequent figures, the bold solid curve 
represents the non-rotating or static limit, and the bold dashed curve the cen- 
trifugal mass shed limit. The long dashed curve is the constant-0 sequence cor- 
responding to the period P = 1.558 ms. The thin solid curves that are roughly 
horizontal are the constant rest mass evolutionary sequences. The evolutionary 
sequences above the maximum stable non-rotating mass configuration, are the 
supramassive evolutionary sequences and those that are below this limit, are 
the normal evolutionary sequences. The almost vertical thin dashed line is the 
limit for instability against quasi-radial modes. The supramassive evolution- 
ary sequences beyond the quasi-radial mode instability limit are represented by 



dotted lines. The numbers against some of the curves in this and all subsequent 
figures represent the baryonic mass for the corresponding sequence; the equilib- 
rium sequence that has 1.4 M & configuration at the static limit is represented 
by an asterisk against it. 

FIG 4: Gravitational mass (Ma) in solar units versus equatorial radius (R) in km 
for EOS model BBB1. 

FIG 5: Neutron star rotation rate (0) in units of 10 4 rad s _1 versus its specific 
angular momentum (cJ/GMq) for EOS model BBB1. The inset shows a close 
up view of the region surrounding the instability limit to quasi-radial mode near 
the mass shed limit. 

FIG 6: Mq - e c relationship for EOS models BBB2. 

FIG 7: M G - R variation for BBB2. 

FIG 8: - cJ/GM* variation for BBB2. 



TABLE 1 



EOS interaction n E /A K E sym (n ) 

BBB1 Avl4+TBF 0.178 -16.46 253 32.5 

BBB2 Paris+TBF 0.176 -16.01 281 32.9 

BPAL21 Skyrme-like 0.160 -16.00 180 30.0 

BPAL32 Skyrme-like 0.160 -16.00 240 30.0 



emp. sat. prop. 0.17 ± 0.01 



-16 ± 1 220 ± 30 30 ±2 



TABLE 2 



EOS e c M G R M M p r orb I 

(g cm~ 3 ) (M ) (km) (M ) (M ) (km) (g cm 2 ) 



BBB1 3.09E+15 

BBB2 3.12E+15 

BPAL21 3.51E+15 

BPAL32 2.67E+15 



1.788 9.646 2.082 

1.917 9.519 2.261 

1.684 9.292 1.940 

1.947 10.509 2.263 



2.356 15.845 1.428E+45 

2.608 16.984 1.593E+45 

2.222 14.921 1.191E+45 

2.579 17.254 1.826E+45 



TABLE 3 



EOS <l c Q. I M G T/W R e u c /tt c J r orb Z p Z f Z b M M v 

BBB1 2.56E+15 1.095 2.428 2.135 0.120 13.129 0.703 0.764 2.658 13.490 0.690 -0.330 1.975 2.471 2.734 

BBB2 2.82E+15 1.203 2.539 2.272 0.123 12.519 0.687 0.825 3.055 13.550 0.849 -0.349 2.483 2.653 3.008 

BPAL21 3.03E+15 1.115 1.904 1.966 0.105 12.604 0.697 0.764 2.123 12.674 0.641 -0.323 1.811 2.253 2.530 

BPAL32 2.27E+15 1.001 3.005 2.300 0.113 14.276 0.699 0.771 3.008 14.611 0.679 -0.328 1.933 2.657 2.962 



TABLE 4 



EOS <l c Q. I M G T/W R e u c /tt c J r orb Z p Z f Z b M M v 

BBB1 2.44E+15 1.079 2.465 2.133 0.120 13.264 0.706 0.756 2.660 13.558 0.677 -0.328 1.935 2.468 2.721 

BBB2 2.82E+15 1.203 2.539 2.272 0.123 12.519 0.687 0.825 3.055 13.550 0.849 -0.349 2.483 2.653 3.008 

BPAL21 2.91E+15 1.100 1.933 1.965 0.106 12.730 0.699 0.757 2.125 12.733 0.630 -0.321 1.779 2.252 2.519 

BPAL32 2.14E+15 0.981 3.070 2.299 0.114 14.481 0.702 0.760 3.011 14.713 0.661 -0.326 1.881 2.655 2.944 



TABLE 5 
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TABLE 5 (contd.) 
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